International Journal of Mechanical and Production 
Engineering Research and Development (IJMPERD) 

ISSN(P): 2249-6890; ISSN(E): 2249-8001 
Vol. 9, Issue 4, Aug 2019, 677-684 
© TJPRC Pvt. Ltd. 

ANALYTICAL AND NUMERICAL APPROACH TO VIBRATION ANALYSIS 
OF LIQUID FILLED ALUMINUM THINCYLINDRICAL SHELLS 

P. KIRAN KUMAR 1 & J. V. SUBRAHMANYAM 2 

'Assistant Professor, Department of Mechanical Engineering, CBIT, Hyderabad, India 
2 Rtd.Professor, Department of Mechanical Engineering, Osmania College of Engineering, Hyderabad, India 

ABSTRACT 

Vibration behavior of empty and variable fluid levels, contained in aluminum cylindrical shells of 2mm and 
3mm subjected to horizontal accelerations are considered. Mathematical expressions, which show the motion of 
cylinder were developed and modified by use of small amplitude wave approximations, enabling equations for the 
various modes ofvibrations and natural frequencies to be obtained [11]. The expression for frequency is formulated by 
the consideration of fluid in the cylindrical shell. The results were compared by considering modal analysis of 
aluminum shells, modeled and analyzed using ANSYS. Natural frequencies for different mode shapes are developed for 
both 2mm and 3mm aluminum shells of empty and with variable water column from the base of the shell. 
Damping ratios were calculated using half power method. Natural frequencies predicted from analytical method 
correlated with ANSYS, were found to be in close agreement. 
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l. INTRODUCTION 

In the recent times, the behavior of storage tanks that carry liquid has got special attention because of 
vulnerability effects of storage tanks due to earthquakes, as they caused liquid spill out due to structural failure, 
environmental effects and leakages, all these are potential losses to the nation. To reduce these effects due to winds 
and seismic excitations, storage tanks are proposed to have mass dampers. In this regard, number of works was 
carried by several researchers, who studied dynamic nature of storage tanks carrying liquids. Most of the works 
fluid mass is added to the structure of the tanks, the dynamic behavior characteristics are obtained with analytical 
approaches by considering simplified geometries and boundary conditions. Some of the earlier works are as 
follows. 

Kubenko and Koval’chuk (2009)[3]observed in their analytical and experimental work that, shells with 
damping models subjected to axial and combined loads, there is big gap in the dynamic behavior, mainly because of 
geometric imperfections. 

Mallon et al 2010[5] used Donnell’s nonlinear shell theory in his analytical work, for the study of 
cylindrical shell with orthotropic material for multimode expansion to discretize PDE to ODE. In this theoretical 
work, shell and electromechanical shaker interaction was considered in the model. They also carried out 
experimental work. There is small gap between experimental and analytical study. 
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Pellicano (2005) [6] conducted experiments when shell of polymeric material with rigid mass on the topis 
subjected to the base excitation. Author claims that, rigid mass undergoes large amplitude while cylinder excited at its first 
mode of vibration with low base excitation about lOg. 

Pellicano, 2007[8]in his further work, hefully considered problem in a linear Point of view. The nonlinear 
phenomena of shell were subjected to base excitation, with sine excitation to the electro mechanical shaker table. Shaking 
the shell due to sine excitation from the bottom induces a vertical motion of the top disk, which causes axial loads due to 
inertia forces. Due to these axial loads, the shell generally subjected to axial-symmetric deformations; but, in some 
conditions, it is noted that a violent resonant phenomenon takes place experimentally, by a strong energy transfer from low 
to high frequencies and large amplitude of vibration. Also, an interesting saturation phenomenon is noticed that the 
response of the top disk was completely flat, as the excitation frequency was changed around the shell base. 

Vijayarachavan and Evan-Iwanowski, 1967 [9] studied parametric instability analysis of circular cylindrical shell 
subjected to seismic excitation, by experimentally and analytically. They considered vertical circular cylinder when the 
base was axially excited, using an electro mechanical shaker. In this study, the in-plane inertia was variable along the shell 
axis, when the base is harmonically excited; it gives rise to a parametric excitation. Experimental and analytical instability 
regions were compared for first axisymmetric mode resonance. In experimental work, glass fiber-reinforced plastic 
cylindrical shell, both empty and water filled were subjected to the nonlinear dynamic deformation of the elastic wall, 
because of hitting by kinematic, two frequency load are studied. 

The present work considers the general problem of thin circular cylindrical shell of 2mm and 3mm subjected to 
horizontal acceleration, containing different levels of in viscid fluid. In order to obtain equations for free surface 
displacements, Small amplitude wave approximations were considered. The fundamental natural frequencies of vibration 
of thin circular cylinder made of aluminum are calculated by determining the added mass of the fluid to the shell structure. 
Different modes of vibration of the cylindrical shell are obtained, with this, added mass of the liquid is added up with the 
structural mass. The natural frequencies were obtained analytically and using finite element approach, and are compared. 

2. MATHEMATICAL FORMULATION 

A. Mathematical Formulation for Determining the Fundamental Natural 
Frequencies of the Circular Cylindrical Shell 

In the view of structural point, while considering the shell expressions, the axial component of motion is not 
considered, and so are the variations of displacement field, axially. Therefore, radial and tangential Love equations are 
written as [1]: 
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Assuming hoop strain to be negligible, Love eqn. in radial direction is reduced to: 
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The mode shapes are of the following admissible type [1]: 

|m w (0) =a M cos(«-0) + p w sin(«-0);v w (0) =a w cos(«-e) +b n sin(«-e)| n= 1,2,. 
This can be conveniently split into two orthogonal families of mode shapes: 

u J ( 0 ) = cos/2-0; vj(d) =-—sin nQ; u^{Q) = sin/2-0; v^( 0 ) =-—cos n 0 \ 

The corresponding mass and stiffness coefficients per unit length are given by 
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The natural frequencies of the structure are given by [1]: 
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Free and rigid in plane translations corresponds to n=l. The shapes constitute a subspace spanned by the two 
orthonormal vectors: 


2 

Ujl = cos0; Vjl =-sin0; u^l^smd^vj =cos0; (9) 

The mode shapes n>l corresponds to axial bending 

B. Considering the Effect of Fluid in the cylinder 

To determine the added fluid mass coefficient inside the cylinder, a small strip model is used. It is assumed in the 
strip model, by considering a narrow strip between z and z+dz, located sufficiently away from the ends z=0 and z= H, we 
now turn our attention towards motion of the fluid forced by a radial vibration of the shell of the type. 


t/(0 : t) = q 1 n cosnQ + q' n sinn6 (10) 

72 = 1 

Where, eql and eq2 stand for modal displacements of the shell Pressure is governed by the boundary value 
problem [1]: 

&E. + ±^E. + J—9+ = 0 (11) 

+ 38 2 
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Splitting the above eqns. into sine and cosine families: 
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Solving the above equations, the pressure field can be determined as: 
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In the present problem, the added mass of the fluid to matrix, as obtained in the vacuum for structural mode and 
the mode shapes of the shell are the same. The common force exerted by the fluid on a shell strip with unit length is 
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Where u. is the unit normal vector in the radial direction, the added mass coefficients per unit length of the shell 
is given by: 


M 




■HR 1 


(19) 


Where, the mode shape is normalized by the condition max u n (1 ’ 2) (0) =1 and m f = pnR 2 H, stands for the physical 
mass of the fluid contained in the shell of unit length. The natural frequencies obtained through this analytical approach 
shown in results and discussion column. 


3. ANSYS FORMULATION 


In modeling section, cylinder is modeled as shell element both 2mm and 3 mm thickness. Finite element model is 
shown in figure-1 for pre processing model. Boundary condition for cylinders for all models fixed at ground position, and 
as well as cap modeled as shell element, as shown in figure-1. Inside, fiuid modeled as fluid element and for fluid structure, 
interaction fluid velocity and density are given for modal analysis and random analysis. For the shell element for the 
Aluminum materials, Structural Solid with 3-D 8-Node - SOLID185 Element and fluid inside the shells as 3-D Acoustic 
Fluid -FLUID30 Element has been considered. Modeled shell with varying water levels (20% level variation from base to 
the shell height) as shown fig-1 has been meshed, using mesh attributes. The dimension of Aluminum hollow shell is of 
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2mm and 3mm thick, lOOmm outer diameter, 300mm height is modeled with lid. The material properties of Aluminum 
with density (p) =2700kg/m 3, modulus of elasticity E =0.683xl0 n N/m 2 , Poisons ratio (0) =0.34, element type SOLID185 
Element is considered for shell element. Fluid with properties of Water density (p) =1000kg/m 3 ' sonic velocity=1498m/s 2, 
Viscosity= 0.00089 pa.s and element type for water as Fluid30, has been considered to prepare model, using ANSYS17.0 

The physical problem has been converted into FEM problem using mapped mesh, and the shell element coupled 
with water element using coupled elements. The geometry is discretized into finite elements using solid 185 and fluid30 for 
structural and fluid, respectively. Outer peripheral nodes of fluid domain and inner peripheral nodes of structural container 
are normally coupled. 



4. RESULTS AND DISCUSSIONS 


The mode shapes for above model, with results are shown in below figure-2&3 for empty and variable water 
column in the shell of 2mm thickness. Modal analysis performed for without water inside cylinder and with water column 
of60, 120,180,240 mm length, from bottom of cylinder. Extracted first four mode shapes and their natural frequencies in 
lanconze method are shown in below figure. Similarly, natural frequencies through modal analysis for 3mm aluminum 
shell were also developed, but mode shapes are not shown with figures but their natural frequencies with damping ratios 
are shown in table 1 &2 for first mode and second mode. It has been observed that acceleration ratio is reduced, with 
increase in thickness, the natural frequencies of the aluminum empty cylindrical shell increases, when water column in the 
cylindrical shell increases, the natural frequencies decrease because of mass of the water level added to the shell. 
And, damping ratios also increase with increases in water column. 
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Figure 2: Four Mode Shapes of Empty Shell 
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Figure 3: (a) 60 mm Water Column Shell Mode Shapes, (b) 120 mm Water Column Shell Mode Shapes, 
(c) 180 mm Water Column Shell Mode Shapes, (d) 240 mm Water Column Shell Mode Shapes 


Table 1: Damping Ratios and Fundamental Natural Frequencies of Aluminum Circular 
Cylinders of 2mm and 3mm Thicknesses for Variable Water Column Heights 


Thickness 

Parameter 

Water Column Length, Mm 

0 

60 

120 

180 

240 

2 mm 

l st Natural 

Frequency Hz 
(Analytical) 

315.4 

276.2 

250.3 

221.8 

178.6 

l s, Natural 

Frequency Hz (FEM) 

318.8 

277.3 

255.3 

215.9 

170.3 

% of error 

1.07 

0.4 

1.99 

2.66 

4.64 

l st mode 

DampingRatio % 

1.14 

2.13 

2.41 

2.60 

2.75 

3 mm 

l st Natural 

Frequency Hz 
(Analytical) 

372.2 

368.1 

351.6 

281.3 

219.0 

l s, Natural 

Frequency Hz (FEM) 

378.4 

375.6 

348.3 

279.4 

224.0 

% of error 

1.66 

2.03 

0.93 

0.67 

2.28 

l st mode 

DampinaRatio % 

2.94 

2.63 

2.81 

2.73 

2.54 


Table 2: Damping Ratios and Second Natural Frequencies of Aluminum Circular Cylinders of 
Thicknesses for 2mm and 3mm Variable Water Column Heights 


Thickness 

Parameter 

Water Column Length, Mm 

0 

60 

120 

180 

240 

2 mm 

2 nd Natural 

F requencyHz( Analy tical) 

431.6 

419.2 

411.4 

400.2 

402.1 

2 nd Natural 

FrequencyHz (FEM) 

438.8 

426.6 

429.6 

410.5 

408.7 

% of error 

1.66 

1.76 

4.42 

2.57 

1.64 

2 nd mode 

DampinaRatio % 

1.50 

3.15 

3.20 

3.38 

3.45 

3 mm 

2 nd Natural FrequencyHz(Analytical) 

501.2 

569.3 

570.1 

538.8 

519.9 

2 nd Natural 

FrequencyHz (FEM) 

504.6 

578.8 

607.5 

550.1 

497.7 

% of error 

0.67 

1.66 

6.56 

2.09 

4.27 

2 nd mode Damping 

Ratio % 

1.71 

3.52 

3.80 

3.95 

4.15 
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CONCLUSIONS 

From table 1&2, it is clear that the natural frequencies of the aluminum empty cylinders increased as thickness is 
increased. The natural frequencies of the empty shell is decreased as the height of the water column increases, due to the 
added mass effect of the water column and presence of the water column in the hollow cylindrical shell. There is increase 
in damping ratios, as length of water column increased. The damping ratios are more for 3mm shell over 2mm thick shell. 
The natural frequencies of the cylinder were determined analytically, and compared with frequencies obtained using 
ANSYS, and are found to be in close agreement. 
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